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ABSTRACT 


The  stress  intensity  factor  K,  a  parameter  used  in  Fracture 
Mechanics,  is  applied  to  the  study  of  the  crack  propagation  stage 
of  fatigue. 

It  was  necessary  to  find  K  for  a  single-edge-notch  specimen 
using  a  fatigue  crack  for  the  notch.  Long  single  wire  strain 
gauges  were  used  to  measure  the  SEN  specimen  compliance  from  which 
K  was  calculated.  Corrections  for  plastic  zone  at  the  crack  tip 
were  applied  to  crack  length.  The  expressions  and  values  for  K 
were  found  to  be  in  good  agreement  with  those  that  had  been 
published  for  use  in  fracture  toughness  tests. 

The  crack  propagation  rate  da/dN,  is  correlated  with  K 
using  a  K  versus  log  (da/dN)  plot  and  a  log  (K)  versus  log  (da/dN) 
plot.  These  correlations  were  affected  by  the  stress  ratio,  and 
the  specimen  thickness.  A  comparison  is  made  with  published 
correlations  and  the  differences  are  discussed.  Fatigue  cracks 
are  discussed  in  relation  to  the  experimental  method  required  to 
measure  the  crack  length,  the  slanting  of  the  fatigue  surfaces, 
the  validity  of  K  for  fatigue  cracks,  and  published  observations  of 
electron  fractography  of  the  fatigue  crack  surfaces. 
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CHAPTER  I 


INTRODUCTION 


1.1  A  Fracture  Mechanics  Approach  to  Fatigue 

Fatigue,  which  may  be  considered  fracture  under  cyclic  loading, 
has  been  one  of  the  most  common  causes  of  service  failures.  Fatigue 
cracks  form  fairly  early  in  the  life  of  materials  after  which  they 
propagate  in  an  assumed  continuous  process. 

The  study  of  the  crack  propagation  stage  of  fatigue  can  be 
aided  by  certain  parameters  of  Fracture  Mechanics.  These  parameters 
are  the  stress  intensity  factor,  which  is  used  in  describing  the 
stress  field  around  the  crack  tip,  and  the  strain  energy  release 
rate  which  is  used  in  describing  energy  changes  of  a  specimen  with 
a  propagating  crack.  Application  of  these  parameters  to  fatigue 
crack  propagation  is  discussed  with  respect  to  corrections  and  the 
use  of  the  single-edge-notch  (SEN)  specimen. 

The  SEN  specimen  has  recently  been  developed  for  studying 
fracture  under  "static"  loadings.  Expressions  of  the  stress  intensity 
factor  were  developed  for  these  fracture  tests.  Studies  on  fatigue 
crack  propagation  have  used  the  SEN  specimen  with  the  expressions 
of  the  stress  intensity  factors  which  were  developed  for  the  crack 
under  plane  strain  conditions.  Experimental  determination  of  the 
stress  intensity  factor  for  the  SEN  specimen  using  a  fatigue 
crack  for  the  notch,  was  thought  useful  in  checking  the  values  of 
the  stress  intensity  factor  of  the  fracture  test.  The  dependence 


of  the  crack  propagation  rate  on  the  stress  intensity  factor  is 
investigated. 
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1.2  Historical  Review 

In  the  past  century,  advances  in  transportation,  automatic 
production  equipment,  and  prime  movers  have  created  a  situation 
where  fatigue  failures  have  become  common.  To  prevent  fatigue 
failures,  the  performance  of  materials  subjected  to  cyclic  loads 
and  stress  concentrations  have  been  investigated  since  the  mid 
1800's. 

In  1870,  Wohler  published  the  results  of  his  twenty  year  study 
on  fatigue  failures.  His  work  resulted  in  the  S-N  curve  which  has 
become  one  of  the  factors  in  evaluating  the  performance  of  eng¬ 
ineering  materials  subject  to  cyclic  loads.  The  use  of  the  S-N 
curve  was  the  basis  for  the  first  scientific  design  of  structures 
subjected  to  repeated  stresses.  This  is  one  of  many  approaches 
that  have  been  used  to  describe  the  basic  mechanism  and  controlling 
factors  of  fatigue,  but  none  have  been  applied  satisfactorily  to 
structural  design, 

Gohn  [1]  and  Hardrath  [2]  have  reviewed  various  fatigue 

theories.  Gohn  states  that 

up  to  a  certain  stage  in  development  each  of  the 
theories  is  in  accord  with  the  observed  phenomena, 
but  to  complete  the  theory  certain  assumptions  are 
usually  made  upon  which  no  experimental  evidence 
is  available.  The  theory,  therefore,  can  not  be 
proved,  but  neither  can  it  be  disproved.! 


1  Gohn,  Fatigue  of  Metals ,  1963,  pp.  107. 
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Considering  the  crack  propagation  stage  of  fatigue,  Paris  and 
Erdogan  [18]  proposed  the  theory  that  the  crack  propagation  should 
be  a  function  of  the  stress  intensity  factor.  They  showed  that  their 
theory  agreed  with  a  wide  range  of  data;  better  than  the  other 
theories  on  crack  propagation. 

The  development  of  the  concept  of  "Fracture  Mechanics"  started 
in  1920,  with  Griffith's  [3]  concept  of  unstable  crack  propagation. 

He  considered  the  crack  length  would  become  unstable  if  the  total 
energy  of  the  system  was  lowered  for  an  increase  in  crack  length. 
Essentially,  he  considered  a  balance  between  the  strain  energy  of 
the  system  and  the  crack  surface  energy.  He  could  then  predict  the 
conditions  under  which  a  flaw  of  given  size  in  the  system,  would  in¬ 
crease  in  size,  usually  resulting  in  complete  failure.  Due  to  evi¬ 
dence  of  plastic  deformation  on  fractured  surfaces,  the  energy 
needed  for  plastic  deformation  was  added  to  the  surface  energy  term 
[4]. 

In  1954,  Irwin  and  Kies  [5]  developed  expressions  for  the  strain 
energy  release  rate  in  terms  of  specimen  compliance  changes  due  to 
crack  growth.  Later,  in  1957,  Irwin  [6]  made  a  theoretical  analysis 
of  the  stress  and  strain  near  the  end  of  a  sharp  crack,  with  the 
stress  intensity  factor  being  the  governing  parameter  of  the  stress 
strain  equations.  The  stress  intensity  factor,  K,  takes  care  of 
the  influence  of  the  configurations,  loads,  and  changing  crack  length 
of  the  specimen  for  the  stress  strain  expressions.  Irwin  [6]  related 
K  to  the  strain  energy  release  rate,  G.  The  stress  analysis  of 
several  specimens  were  given,  from  which  K  was  calculated. 
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For  specimens  of  an  arbitrary  design,  alternative  methods  will 
yield  sufficiently  accurate  K  expressions.  These  are: 

(1)  Theoretical  stress  analysis  procedures  [13,  14] 

(2)  Experimental  strain  analysis  procedures  using  photoelasticity 
and  strain  gauges  [6] 

(3)  Experimental  measurements  of  specimen  compliance  as  a  function 
of  crack  length  [5], 

Since  World  War  II  the  extensive  investigations  of  many  fractures 
has  led  to  most  of  the  development  and  application  of  Fracture  Mech¬ 
anics  to  the  problems  associated  with  fatigue  crack  propagation. 

R.E.  Peterson  [7],  in  reviewing  design  aspects  required  in  preventing 
fatigue  failures,  cites  the  widely  known  fatigue  failures  of  the 
Comet  I  aircraft  as.  a  major  incentive  in  the  aircraft  industry  to 
investigate  fatigue  and  design  to  prevent  such  failures  in  larger 
modern  aircraft.  Consequently,  most  of  the  recent  work  in  fatigue 
crack  propagation  has  been  done  by  research  facilities  affiliated 
with  the  aircraft  industry. 

In  the  early  nineteen  fifties,  a  steel  with  a  yield  strength 
above  200,000  psi  had  been  developed  and  was  selected  for  use  in  a 
Polaris  Motor  Case  [8].  Some  of  these  missiles  unexpectedly  blew 
up  in  the  air,  others  burst  in  hydrostatic  tests  and  some  even 
fractured  in  manufacturing  due  to  the  brittle  nature  of  the  high 
strength  steel.  As  a  result,  in  1958  the  special  ASTM  Committee  on 
Fracture  Testing  of  High-Strength  Sheet  Material  was  set  up  to  pro¬ 
vide  a  focal  point  of  interest,  talent,  and  information  in  the 
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field  of  fracture..  The  wide  acceptance  and  use  of  K  and  G  was  partly 
due  to  the  five  reports  of  the  committee. 

The  understanding  of  the  strength  and  fracture  of  materials  has 
increased  tremendously  in  the  past  decade  both  from  the  macroscopic 
(fracture  mechanics)  viewpoint  and  the  microscopic  viewpoint.  Con¬ 
ferences,  such  as  the  International  Conference  of  Fracture  [9]  held 
in  Sendai,  Japan  in  September  12-17,  1965  which  presented  both  the 
macroscopic  and  microscopic  viewpoint  tended  to  promote  more  complete 
understanding  of  fractures.  The  concepts  of  fracture  of  solids  has 
advanced  to  the  point  where  they  can  be  used  by  practising  engineers. 

With  the  fairly  recent  development  of  the  fracture  parameters 
K  and  G,  many  test  methods  and  specimens  have  been  used  to  measure  K 
and  G  for  various  materials.  The  SEN  specimen  has  been  one  of  the 
more  recent  specimens  developed  to  enable  testing  with  smaller 
machines  and  less  equipment. 

1.3  History  of  the  Single-Edge-Notch  (SEN)  Specimen 

The  SEN  specimen  was  initially  introduced  by  Sullivan  [11]  as 
a  new  specimen  design  for  plane  strain  fracture  toughness  tests. 

It  has  definite  advantages  over  the  center  notch  and  double  edge 
notch  sheet  specimens  that  were  currently  being  used.  Less  material, 
smaller  loading  capacity  of  the  test  machines  are  required,  and 
only  one  crack  tip  is  observed  in  comparison  with  two  crack  tips 
of  other  specimens.  The  SEN  specimen  is  a  metal  sheet,  containing 
a  crack  (notch)  at  the  edge  midway  between  the  supports,  which 
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propagates  from  one  edge  of  the  sheet  to  the  other;  due  to  a  pin 
loading  applied  at  the  ends  of  the  long  specimen.  The  continuously 
changing  configuration  and  stress  distribution  of  the  SEN  specimen, 
brought  about  by  the  growing  crack,  are  accounted  for  by  the  stress 
intensity  factor. 

Sullivan  [11]  employed  experimental  compliance  measurements 
to  obtain  expressions  for  K.  The  small  size  of  the  test  specimen 
and  certain  features  of  her  test  procedure  caused  Srawley,  Jones,  and 
Gross  [15]  to  question  the  accuracy  of  the  K  values.  They  found  more 
accurate  values  of  K  using  larger,  longer  specimens  together  with 
more  accurate  experimental  techniques.  In  order  to  check  the  accur¬ 
acy  of  these  values,  Gross,  Srawley,  and  Brown  [13]  applied  boundary 
collocation  procedures  to  a  suitable  stress  function,  to  obtain  a 
K  expression  for  the  SEN  specimen.  Srawley,  Jones,  and  Gross  [15] 
found  their  K  expressions  to  be  in  excellent  agreement  with  those 
found  by  Gross,  Srawley,  and  Brown  [13],  and  disagreeing  with  the 
K  expressions  of  Sullivan's  [11]. 

The  K  expressions  of  references  13  and  15  were  for  a  sharp 
crack  in  a  perfectly  elastic  material  in  a  plane  strain  state. 

However  an  actual  crack  has  a  plastic  zone  at  the  crack  tip  for  which 
a  correction  must  be  made. 

The  SEN  specimen  is  now  commonly  being  used  in  fatigue  crack 
propagation  studies.  In  this  investigation  it  was  thought  useful 
to  obtain  K  values  using  a  SEN  specimen  with  an  actual  fatigue  crack 
for  the  notch,  rather  than  an  artificial  crack  as  used  by  Srawley 


et  al . 


CHAPTER  II 


THEORY 


2.1  Griffith's  Work 

Griffith  [3]  proposed  a  crack  propagation  theory  in  1920 
that  was  the  ground  work  of  the  present  fracture  mechanics  energy 
concepts.  He  investigated  the  effects  of  surface  treatments  on 
the  strength  of  metallic  machine  parts.  Griffith  found  that  small 
scratchs  reduced  the  strength  of  parts  and  that  the  hypothesis  of 
rupture  of  the  present  day  did  not  account  for  this  effect. 

In  considering  a  theoretical  criteria  of  rupture  he  adopted 

the  Theorem  of  Minimum  Potential  Energy,  and  applied  it  to  an 

elliptic  shaped  flaw  in  a  plate.  His  statement  of  the  criterion  is, 

if  a  system  can  pass  from  the  unbroken  to  the  broken 
condition  it  must  do  so  by  a  process  involving  a 
continuous  decrease  in  potential  energy. 2 

If  the  potential  energy  of  a  specimen  containing  a  crack  is  differ¬ 
entiated  with  respect  to  the  crack  length,  and  if  the  differential 
is  positive,  there  is  an  increase  in  potential  energy  for  an  increase 
in  crack  length.  The  crack  is  therefore  stable.  If  on  the  other- 
hand  the  differential  is  negative,  the  potential  energy  would  de¬ 
crease  for  an  increase  in  crack  length  and  hence  the  crack  is  un¬ 
stable  and  will  increase  in  length  until  the  potential  energy  has 
reached  a  minimum  value. 


2  Griffith,  Ruptures  and  Flow  in  Solid,  Philosophical  Transaction, 
Vol .  221,  p.  165. 
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Griffith  considered  the  changes  in  the  potential  energy  for  a 
specimen  with  a  crack  to  be  made  up  of  the  increase  in  strain  energy 
W,  and  the  surface  energy  S.  To  illustrate  Griffith's  concept, 
consider  an  example  from  Weiss  and  Yukawa  [4]  where  an  elastic  crack 
of  length  2a  in  an  infinite  plate  of  unit  thickness  is  considered. 

A  uniaxial  stress  a  is  applied  at  infinity,  perpendicular  to  the 
crack.  The  total  change  of  strain  energy  due  to  the  crack  is 

,,  -o2naz 

w - — 

The  energy  gain  of  the  plate  due  to  the  creation  of  the  new  surface 
having  a  surface  tension  T,  is 

S  =  4aT 

The  total  decrease  of  potential  energy  is  W  plus  S.  If  the  crack 
length  is  taken  as  the  only  variable,  the  conditions  for  equilibrium 
and  stability  are  found  by  differentiating  the  potential  energy  with 
respect  to  the  crack  length.  Let  the  initial  value  of  potential 
energy  for  zero  crack  length  equal  VQ. 

The  differential  is 

^(w  +  s  +  v0)  (1) 

If  then  ^  (W  +  S  +  V  )  <  0  (la)  The  crack  length  is  stable 

=  0  (lb)  The  crack  length  is  critical 

>  0  (lc)  The  crack  length  is  unstable 

Substituting  the  values  of  W  and  S  in  equation  (lb)  yields 
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as  the  value  of  stress  at  which  the  crack  length  is  on  the  verge 
of  increasing  in  length  and  is  considered  the  critical  crack  length. 

Experimental  work  of  the  nineteen  forties  showed  plastic  de¬ 
formation  to  be  evident  on  fracture  surfaces.  Orowan  concluded  that 
the  changes  in  potential  energy  due  to  crack  growth  must  include  the 
work  done  in  plastic  deformation  [4].  The  absolute  change  in  the 
strain  energy  must  equal  the  absolute  changes  in  surface  energy  plus 
the  plastic  work.  The  change  in  strain  energy  for  an  increment  in 
crack  growth  for  a  unit  thickness,  is  known  as  the  strain  energy 
release  rate  G.  The  strain  energy  release  rate  can  be  calculated 
by  theoretical  or  experimental  analysis.  The  other  parts  that  make 
up  the  changes  in  strain  energy  are  unknown.  Values  of  G,  found 
experimentally,  can  be  used  to  predict  the  conditions  for  which  a 
crack  length  will  become  unstable.  These  are  known  as  the  critical 
values  of  G  or  the  fracture  toughness  G  .  For  the  example  considered 

V 

previously,  the  expression  for  the  energy  release  rate  is: 

G  =  (2) 

In  the  following  sections,  G  is  related  to  the  stress  strain 
field  around  a  crack  tip  and  to  experimentally  measurable  values  of 
the  compliance  of  specimens.  These  relationships  are  used  for  the  SEN 
specimen  under  consideration. 
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2.2  Stress  Strain  Field  Around  the  Crack  Tip 

Irwin  [6]  obtained  expressions  for  the  stress  strain  field 
around  the  crack  tip.  He  used  a  semi-inverse  method  developed  by 
Westergaard  [16]  for  solving  certain  classes  of  plane  strain  and 
plane  stress  problems.  Paris  and  Si h  [14]  illustrate  this  semi¬ 
inverse  method  in  detail. 

Irwin  assumed  the  material  to  be  a  homogeneous,  isotropic, 
elastic  solid  with  an  infinitely  sharp  crack  with  a  straight  crack 
front  perpendicular  to  the  specimen  surface,  (a  flat  crack).  The 
expressions  for  the  stress  field  around  the  crack  tip  shown  in 
Figure  2.1  are  listed  as  follows: 


a 


x 


K  .  0  .  30 \  /  o  \ 

-  cos  2  (1  -  sin  J  sin  (3a) 

v^Trr 


a 


y 


K  0  /  -i  ,  .  0  30  \  /  o.  \ 

-- —  cos  J  (1  +  sin  7,-  sin  -y)  (3b) 

/ZttF 


(3c) 


The  displacements  are: 


u 


v  =  jr  (1  +  y )(■"”)  sin  (2  -  2y  -  cos2  |-)  (4b) 


In  the  development  of  the  above  expressions  the  terms  containing  r/a 
have  been  considered  negligible  and  are  neglected.  Thus  these 
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Figure  2.1  Coordinates  For  the  Region  Around  a  Crack  Tip 


Figure  2.2  Increment  of  Crack  Growth 
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expressions  are  for  sheet  problems  and  apply  only  to  the  region 
around  the  crack  tip, 

2.3  Relating  the  Strain  Energy  Release  Rate  to  the  Stress  Intensity 

Factor 

Irwin  [6]  introduced  a  method  for  calculating  G  from  the  elastic 
stress  analysis  of  a  crack  tip.  Since  the  region  around  a  crack  tip 
is  completely  elastic,  the  decrease  in  strain  energy  due  to  an  incre¬ 
ment  of  crack  growth,  6a,  is  equal  to  the  increment  of  work,  6W, 
required  to  close  the  increment  of  crack  growth.  Consider  an  elastic 
body  of  unit  thickness,  containing  a  crack  of  length  a  which  has 
been  extended  by  6a  as  shown  in  Figure  2.2.  Two  sets  of  axis  y'x  and 
yx  are  needed  to  describe  conditions  for  the  cracks  of  length  a  and 
a  +  a. 

The  loss  of  strain  energy  due  to  the  crack  growth  6a  will  be 

restored  by  pulling  the  crack  surfaces  of  6a  together.  The  crack 

surfaces  are  stress  free  and  the  crack  surface  displacement  v  is 

given  by  the  equation  (4a),  where  the  distance  r  is  related  to  the 

xy  origin.  Consider  the  element  of  surface  area,  1  •  dr,  moving  a 

distance  v  with  the  force  on  this  area,  a  dr,  where  a  increases 

linearly  from  zero  to  a  maximum,  as  the  displacement  is  accomplished. 

The  maximum  value  of  a  is  known  from  the  equation  (3a),  where  the 

«y 

crack  is  of  length  a.  Note  that  the  displacement  v  is  referred  to 
the  xy  axis  and  is  a  function  of  the  radius  r.  The  maximum  value 
of  stress  a  is  referred  to  the  xy'  axis  and  is  a  function  of  the 

<y 
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radius  (a-r). 

The  work  required  to  close  the  crack  (both  surfaces)  over  the 
distance  a  is 

a  v 

(5a) 


o  v 

<SW  =  2  /  “  -i-  dr 

J  o  2 


where,  a  refers  to  the  xy"  axis  with  0=0  and 


K 


^  A  ( a-  r ) 


The  displacement  v  refers  to  the  xy  axis  with  0  =  180°  and  from 
equation  (4a) 

V  =  %h!1  K  ri/z 

TT  E 


(6) 


Thus , 


Ga  =  6W 


(5b) 


6W  =  2  /  a  -Kv  2 ■■■-■  (1  -  v2)  dr 

*  F  tt  ^  ' 


(5c) 


Integrating  yields : 


Ga 


2K'  (1  -  U2 )  r-OlTTn 

=  1 - E - LTJ 


(5d) 


The  relations  of  the  strain  energy  release  rate  G  in  terms  of  the 
stress  intensity  factor  for  the  plane  strain  state  becomes 

.  K2n  -  u2) 

E 


G 


(7a) 


« 
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The  relation  for  the  plane  stress  state  can  be  found  similarly  and 
is 

G  =  Y-  (7b) 

As  a  result  of  the  above  development,  the  strain  energy  release 
rate  G  is  also  known  as  the  crack  extension  force. 

2.4  Relating  the  Strain  Energy  Release  Rate  to  the  Compliance  of 
the  Test  Specimen 

Irwin  and  Kies  [5]  related  the  strain  energy  release  rate  G  to 
the  change  of  compliance  of  the  specimen.  A  more  general  development 
of  the  relationship  is  given  as  follows. 

Consider  the  specimen  shown  in  Figure  2.3  with  an  initial  load 
of  PQ  and  crack  length  a  .  The  displacement  e,  between  the  loading 
points  increases  proportionally  with  the  load  as  shown  in  Figure  2.4. 
Compliance,  C,  is  defined  as  the  displacement  eQ  divided  by  the  load 
P  and  is  actually  the  inverse  of  stiffness. 

During  an  increase  of  crack  length  6a  the  load  would  decrease  an 
amount  6P,  and  the  loading  points  extend  an  amount  6e.  The  total 
strain  energy  before  the  crack  extension  was  P  C  and  after  is 
(P  -  6P)(eo  +  6e).  The  total  change  in  strain  energy,  6V,  is  denoted 
in  Figure  2.4  by  the  horizontal  cross-sectional  area  minus  the  angular 
cross-sectional  area.  The  horizontal  cross-sectional  area  is  equal 
to  the  increment  of  work,  6W,  done  on  the  system,  by  the  external  force 
P.  The  strain  energy  released,  6aG,  due  only  to  the  increment  of  crack 
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Figure  2.4  Total  Strain  Energy 
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growth,  6a,  is  equal  to  the  total  strain  energy  change  minus  the 
work  done  by  external  forces,  as  given  in  the  following  expression: 


G6a  =  6W  -  6V 


G  is  defined  for  an  infinitesimal  element  of  crack  growth  6a.  Thus, 


G  =  1  i  m 

6a  ■+  0 


6W  -  6V 
6a 


(8) 


or 


r  _  dW  dV 
b  da  “  da 


Total  strain  energy  expression  is 


V  = 


Using  the  compliance  relation  (e  =  CP), 


V  = 


CP: 


d V  =  CP3P  IP2  3C 
da  3a  2  3a 


(10) 


The  work  term  is  6W  =  P6e0  Differentiating  with  respect  to  crack 
growth  results  in: 


dW 

3a 


P3e 
3  a 


PC3P  +  P23C 


3a 


3  a 


(ID 


Substituting  equation  (10)  and  (11)  in  equation  (8)  yields 


G  = 


P2  3C 

IT  3a 


(12) 
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Equation  (12)  indicates  that  G  is  independent  of  small  changes  in 
load  or  loading  points,  and  is  dependent  only  on  the  changes  of 
specimen  compliance  for  changes  in  crack  length  and  load.  Experi¬ 
mental  values  of  compliance  can  be  used  in  equation  12  to  determine 
expression  for  G.  The  compliance,  measured  at  different  crack  lengths 
can  be  plotted  on  a  graph  against  the  crack  lengths.  The  slope  of 
this  plot  is  used  to  determine  G, 

2.5  Corrections  for  Plasticity  Effects 

A  small  zone  of  yielding  is  always  present  at  the  crack  tip. 

If  the  plastic  zone  is  "small"  and  the  nominal  stress  in  the  specimen 
is  less  than  eighty  percent  of  the  yield  stress,  the  elastic  analysis 
may  be  applied  with  small  corrections  for  the  crack  length  [17],  In 
fracture  toughness  studies,  it  has  been  acceptable  to  make  corrections 
for  the  small  plastic  zone.  These  corrections  also  should  be  applied 
in  fatigue  work. 

Consider  a  plastic  zone  model  shown  in  Figure  2.5  where  the 
yielded  region  around  the  crack  tip  tends  to  blunt  the  crack  tip. 

The  crack  surfaces  would  move  apart  with  no  further  crack  growth,  as 
suggested  in  Figure  2.5.  The  elastic  region  around  the  crack  tip 
would  "see"  an  elastic  sharp  crack  tip  at  a  different  point  than  the 
yielded  blunted  crack  tip.  For  a  blunted  crack  length  of  a,  the 
elastic  region  around  the  yielded  zone  would  "see"  an  elastic  crack 
length  of  a  +  r  .  The  elastic  stress  field  in  the  region,  a  little 
distance  away  from  the  actual  blunted  crack  tip  of  length  a  is  the 
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Figure  2C5  Plastic  Zone  Model 
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same  as  that  for  the  ideal  elastic  crack  of  length  a  +  r  . 

In  the  elastic  analysis,  the  yield  stress  would  occur  at  a 
radius  r  ,  calculated  from  equation  (3b)  for  0  equal  to  zero  and  a 

y  ] 

equal  to  the  yield  stress,  a  ,  where 

y* 


r  =  _L  (JL)2 

y  2tt  QyS 


(13) 


If  the  yielded  zone  is  considered  to  be  a  cylinder  as  shown  in 
Figure  2,5,  then  the  apparent  elastic  crack  length  would  be  at  the 
center  of  the  yielded  zone,  making  the  apparent  elastic  crack  length 
a  +  r  .  The  actual  shape  of  the  yielded  zone  is  not  known,  but  this 
model  is  an  acceptable  approximation  [17]„ 

The  yield  stress  for  aluminum  alloys  is  usually  taken  as  the 
.2%  strain  offset  of  the  stress  strain  curve.  This  yield  stress  is 
taken  from  a  tension  specimen  in  a  state  of  plane  stress,  which  would 
be  lower  than  the  yield  stress  taken  from  a  specimen  in  a  state  of 
plane  strain.  Irwin  [17]  indicated  that  the  plastic  zone  along  a 
plane  strain  segment  of  the  crack  border  may  be  smaller  than  that  of 
the  plane  stress  portion  because  of  the  change  in  shear  stress  and 
other  considerations.  He  applied  a  reduction  factor  of  3  to  r^  of 
equation  (13).  Hence  for  plane  strain  problems,  he  suggests 


_L  (. 

6tt  a 


K 


ys 


(14) 


Crack  lengths,  measured  by  the  microscope  are  blunted  and  must 
be  corrected  to  equivalent  elastic  cracks  before  applying  elastic 
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analysis .  The  blunted  crack  length  may  initially  be  used  to  obtain 
an  initial  value  of  K„  The  plastic  zone  correction  r  ,  may  be  then 
calculated  using  equations  (13)  or  (14)  and  the  apparent  elastic  crack 
may  be  found.  An  iterative  procedure  may  be  used  to  determine  K. 

2.6  Fatigue  Crack  Propagation 

Fatigue  has  been  defined  as  fracture  under  cyclic  loading  and  is 
the  most  common  cause  of  service  failures  [10],  A  fatigue  failure 
can  be  thought  to  have  three  stages:  initiation,  propagation,  and  a 
final  unstable  failure.  Initiation  is  the  period  required  for  a  crack 
to  start  growing  from  a  micro  crack  or  structural  flaw  to  a  propaga¬ 
ting  crack.  In  the  second  stage  the  crack  grows  slowly  and  almost 
continuously  due  to  the  cyclic  loading.  The  final  stage  is  rupture 
due  to  an  unstable  running  crack  and/or  ductile  yielding. 

Various  approaches  have  been  used  to  explain  the  crack  propaga¬ 
tion  stage,  resulting  in  an  abundance  of  crack  propagation  laws. 

Paris  and  Erdogan  [18]  in  1964  did  a  critical  analysis  of  the  more 
acceptable  crack  propagation  laws  and  their  approach  agreed  better 
than  the  other  crack  propagation  laws  for  a  wide  range  of  test  data. 
Paris  and  Erdogan  proposed  that  the  crack  propagation  rate  should  be 
a  particular  function  of  the  stress  intensity  factor  for  a  particular 
material.  A  plot  of  the  logarithm  ( AK)  versus  the  logarithm  ( da/dN ) 
was  used  to  show  that  there  was  general  agreement  for  the  law, 


m  a(AK)4 


(15) 
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The  more  recent  paper,  Hertzberg  and  Paris  [19]  uses  the  stress  in¬ 
tensity  approach  along  with  the  results  of  electron  fractography 
studies  in  studying  crack  propagation„  They  showed  that  the  striations 
observed  by  the  electron  microscope  had  almost  the  same  correlation 
with  AK,  as  the  crack  propagation  rate  on  log-log  plots.  In  crack 
propagation  studies,  there  is  not  substantial  agreement  of  AK  versus 
da/dN  correlations  from  different  sources.  The  log  (AK)  versus  log 
( da/dN )  plot  minimizes  the  differences  of  the  data  so  a  single  line 
may  be  drawn  through  the  data.  One  reason  for  the  differences  is  that 
many  factors  which  affect  the  crack  propagation  rates  are  not  accounted 
for  by  the  single  parameter  AK.  The  stress  intensity  factor,  AK,  is 
based  only  on  the  value  of  the  cycling  stress  and  specimen  geometry. 
Other  factors  affecting  the  crack  propagation  rate  are  the  stress  ratio 
(R),  type  of  the  waveform  of  the  cyclic  load,  frequency  of  the  cycling 
load,  size  of  specimen  (thickness),  temperature,  and  even  such  things 
as  the  humidity  of  the  test  room. 

Roberts  and  McEvily  [20]  have  plotted  log  (aK)  versus  log  (da/dN) 
for  various  values  of  R.  They  concluded  that  the  crack  propagation 
rate  is  a  function  of  AK  and  R  as  follows: 

^  a  (R)'2  (AK)1*  (16) 

2.7  Application  of  Theory  to  Test  Data 

The  change  in  compliance  due  to  the  crack  was  found  experimentally 
for  25  to  50  crack  lengths  varying  from  0  to  1.7  inches  on  the  3.5  inch 
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wide  SEN  specimen.  The  strain  energy  release  rate  G,  was  found  from 
the  compliance  measurements  using  equation  12: 

r  _  P2  dC 

G  =  T  dl 

where  P  is  the  applied  force  per  unit  thickness,  and  C  is  the  total 
compliance  of  the  specimen  per  unit  thickness.  It  is  convenient  to 
express  this  relationship  in  the  dimensionless  form 

EwG  E  dC  fT-j\ 

p  2  "  2  da/w  [U) 

The  dimensionless  form  of  the  strain  energy  release  rate  as  in 
equation  (17)  may  be  correlated  with  the  relative  crack  length  a/w. 

This  correlation  will  apply  to  specimens  of  any  material  or  size, 
providing  the  length-width  ratio  of  specimens  are  within  certain  limits. 

In  calculating  the  stress  intensity  factor  from  the  strain  energy 
release  rate,  there  is  a  question  of  whether  the  two-dimensional  plane 
stress  or  plane  strain  analysis  may  be  used  for  test  specimens.  The 
state  of  stress  in  a  cracked  plate  specimen  is  one  of  plane  stress, 
except  within  a  very  restricted  region  of  the  crack  where  the  con¬ 
straining  influence  of  the  crack  induces  a  stress  state  approaching 
plane  strain.  The  crack  surface  and  the  region  near  the  center  of 
the  crack  tip  is  considered  to  be  in  plane  strain  state  for  lower  values 
of  K  and  flat  cracks.  The  plane  stress  state  appears  to  start  at  the 
surface  of  the  specimen  with  the  formation  of  a  shear  lip.  The  shear 
lip  gradually  moves  across  the  crack  tip  to  the  center  of  the  specimen 
for  increased  values  of  K.  It  is  generally  agreed  that  once  the  sur¬ 
face  of  the  crack  attains  a  slope  of  45  degrees,  the  region  around  the 
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crack  tip  is  in  a  state  of  plane  stress.  Irwin  [21]  has  discussed 
the  question  of  plane  stress  or  strain  for  fracture  toughness  tests. 
Hertzberg  and  Paris  [19]  use  some  of  Irwin's  arguments  to  discuss 
fatigue.  The  stress  intensity  factor  has  been  calculated  for  both 
plane  stress  and  plane  strain  using  equations  (7a)  and  (7b).  The 
factor  K  has  been  developed  for  a  perfectly  elastic  material  and  since 
a  plastic  zone  is  always  present  at  the  crack  tip,  the  apparent  elastic 
crack  lengths  must  be  calculated. 

In  these  experiments,  a  microscope  was  used  to  measure  the  surface 
crack  lengths.  The  plastic  zone  correction,  r  of  equation  (13)  was 
added  to  the  surface  crack  lengths  for  the  plane  stress  problem.  For 
the  plane  strain  problem,  the  crack  length  at  the  center  of  the  crack 
front  was  obtained  by  adding  the  correction  for  the  parabolic  shape 
of  the  crack  front  to  the  surface  crack  length.  The  correction  was 
found  to  be  zero  at  the  start  of  the  fatigue  crack  and  increased  to 
maximum  value  of  0.005,  0.007,  and  0.009  inches  at  the  end  of  the 
tests  for  specimen  thicknesses  of  0.050,  0.080  and  0.126  inches  respec¬ 
tively.  The  plastic  zone  correction  of  equation  (14)  was  added  to  the 
center  crack  lengths  for  the  plane  strain  problem. 

Since  K  depends  on  the  crack  length  and  the  crack  length  correction 
depends  on  K,  an  iterative  process  was  used  to  calculate  K.  The  stress 
intensity  factor  was  calculated  using  the  original  measured  crack 
lengths.  The  plastic  zone  corrections  were  calculated  using  the  initial 
K  value  and  applied  to  the  crack  lengths.  A  new  K  was  then  calculated 
using  the  corrected  crack  lengths.  The  process  was  repeated  until  it 
was  evident  that  no  further  iterations  were  necessary. 
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The  plot  of  compliance  measurements  versus  relative  crack  length 
were  fitted  to  a  polynomial  expression  of  a/w  of  degrees  3  to  10  by 
means  of  a  least-square-best- fi t  digital  computer  program.  The  result¬ 
ing  expression  was: 


EAC  =  4(i)2+^(£)3+^(£)^  + 


2  /  a \ 3  .  "3 


n 

rTTT 


a\n+l 


(-) 

vw' 


for  n  =  4,  5, 


...  i 


10 


(18) 


The  constant  and  first  power  of  a/w  were  not  included  in  the  expression 
because  the  strain  energy  release  rate  is  zero  for  zero  crack  length. 
The  criteria  for  optimum  fit  [22]  was  taken  to  be  the  sum  of  the  re¬ 
siduals  squared,  divided  by  the  degree  of  freedom  m  -  k  -  1,  where  m 
is  the  number  of  points  and  k  is  the  number  of  terms  in  the  polynomial. 
This  error  variance  indicates  that  polynomials  of  degrees  3  to  5  were 
usually  as  good  as  higher  degree  polynomials  for  the  single  precision 
computer  program. 

Using  equations  (7b),  (12),  (13),  and  (18)  the  stress  intensity 
factor  for  the  plane  stress  was  found  from  compliance  measurements  and 
was  put  in  the  following  form: 

K  =  2^L  [A]  !£+  A,h£)2  +  ...  +  A„(^-)n]l/2  (19) 


w 
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n'  w 


where 


1  ,  K  >2 


a_  =  a  +  Ty-  (-±-) 

2tt  a 

ys 


(20) 


and  stress  intensity  factor  for  the  plane  strain  case  is: 

51  51  51 

K  =  gy/w  r  A  _SL  +  A  +  +  A  ( j£.lnl  l/ 2 

on  2)  LA1  w  VV  +  +  V  w;  J 

2(1 -y  ) 

1  K 

where  a  =  a  +  ^—  ( - )2  +  curvature  correction. 

c  6TT  oys 


(21) 
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The  above  equations  are  presented  graphically  in  Chapter  4. 
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CHAPTER  III 


EXPERIMENTAL  TEST  METHODS 


3.1  Test  Specimen 

The  material  used  in  this  investigation  was  bare  aluminum 
alloy  2024-T3  which  is  one  of  the  common  materials  used  in  crack 
propagation  studies  [2,  19,  20].  Hence  a  comparison  of  results  can 
be  made  quite  easily.  The  clad  material  2024-T3  is  presently  used 
in  large  commercial  quantities,  but  the  bare  2024-T3  was  chosen 
because  of  the  effects  the  soft  cladding  may  have  on  crack  propagati 
results.  The  mechanical  properties  of  the  material  listed  in 
Table  I  were  obtained  from  the  tension  specimen  shown  in  Figure  3.1. 

The  SEN  specimens  shown  in  Figure  3.2  were  cut  from  three 
sheets  of  material  of  thicknesses  0.050,  0.080,  and  0.125  inches. 
They  were  cut  using  a  milling  cutter  with  the  direction  of  rolling 
parallel  to  the  applied  stress  and  perpendicular  to  the  path  of  the 
propagating  crack.  The  length  to  width  ratio  for  SEN  specimens 
should  be  greater  than  four.  Srawley,  Jones,  and  Gross  [15]  found 
no  significant  differences  between  K  curves  obtained  from  specimens 
with  a  length  to  width  ratio  of  four  and  specimens  with  a  length  to 
width  ratio  of  eight.  Compliance  readings  for  these  specimens  are 
shown  in  Figure  4.4  Compliance  readings  on  the  24  inch  specimen 
were  about  three  percent  lower  than  those  of  the  12  inch  specimen. 
This  difference  was  not  considered  to  be  significant. 
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TABLE  I 

Mechanical  Properties  of  Bare  2024-T3 


Mechanical  Properties 
of  Aluminum  Alloy 

Test  Specimen  Thickness 

0.05  Inches  0.080  Inches  .125  Inches 

Gy  ( .2%  Offset)  k.s ,i . 

55.5 

52.0 

57. 

a  (ultimate)  k.s .i . 

71.5 

67.8 

74.8 

E  k .  s .  i . 

10,200 

10,360 

10,300 

y 

.320 

.326 

.325 

%  elongation 

(6  inch  gauge  length 

22.1 

21.2 

11.0 

Figure  3.1  Tension  Test  Specimen 
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Figure  3.2  SEN  Specimen 
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An  ASTM  Special  Committee  [23]  recommended  types  of  holding 
blocks  for  test  specimens,  but  the  limited  length  of  testing  space 
in  the  fatigue  machine  (15  inches)  required  a  special  design  of 
holding  fixtures. 

Two  quarter  inch  steel  plates  were  bolted  to  each  end  of  the 
specimen.  Loads  were  applied  by  two,  half  inch  diameter  bolts  through 
the  center  of  the  steel  plates.  The  plates  were  held  in  the  fatigue 
machine  so  as  to  permit  free  rotation  about  the  bolts,  but  preventing 
the  specimen  from  twisting.  The  effect  of  friction  on  the  rotation 
of  the  specimen  may  have  been  eliminated  by  the  use  of  special 
bearings,  but  this  requirement  was  not  deemed  necessary. 

The  area  of  the  surface  of  the  specimen,  to  which  the  strain 
gauges  were  glued,  was  sanded  lightly  and  cleaned  to  insure  a  good 
bond  of  the  strain  gauges  to  the  surface.  Care  was  taken  not  to 
scratch  or  remove  the  oxide  film  from  the  surface  of  the  specimen 
where  the  crack  was  to  propagate.  The  darker  oxide  film  on  the 
surface  of  the  specimen  made  the  shiny  crack  surface  easy  to  see. 

3.2  Testing  Machine  and  Room 

A  Sonntag  Universal  Fatigue  Testing  Machine,  Model  SF-l-U, 
with  an  alternating  load  frequency  of  1800  cycles  per  minute  was 
used.  Calibration  tests  indicated  a  static  load  error  less  than 
three  percent  of  load  or  one-half  percent  of  the  total  capacity. 

There  were  other  sources  of  errors  in  loads  besides  the  errors 
indicated  in  calibration  tests.  The  static  load  was  set  manually 
at  the  start  of  each  test.  The  static  load  changed  only  slightly 
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during  the  test  because  of  the  compliance  changes  of  the  specimen,  thus 
the  static  load  was  not  reset  during  the  test.  A  deflection  of  the 
head  fixtures  of  the  testing  machine  was  noted  when  the  dynamic  load 
approached  the  capacity  of  the  machine.  The  magnitude  of  the  load 
errors  due  to  the  deflection  of  the  head  fixture  and  the  manual  load 
setting  were  unknown  but  these  should  not  be  significant. 

The  testing  machine  was  located  in  a  room  without  air-conditioning. 
The  tests  were  carried  out  in  June  and  July  with  the  temperature  re¬ 
maining  fairly  constant,  with  a  maximum  of  80  degrees  and  a  minimum  of 
73  degrees  during  testing.  The  relative  humidity  was  found  using  a 
wetbult  thermometer,  and  remained  between  25  to  45  percent  during  all 
tests.  Although  the  temperature  and  humidity  were  not  controlled, 
daily  variations  were  within  two  degrees  of  temperature  and  five  per¬ 
cent  humidity.  These  daily  variations  were  considered  small  and 
would  not  make  any  significant  differences  in  the  tests. 

3.3  Test  Loads 

A  maximum  stress,  based  on  the  gross  cross-sectional  area, 
of  4.5  ksi  and  a  stress  ratio,  R,  equal  to  .2  was  applied  to  all 
specimens.  An  alternating  stress  of  3.6  ksi  was  chosen  so  the 
crack  would  propagate  across  the  specimen  in  a  reasonable  time 
(4  to  6  hours).  The  minimum  stress  was  set  as  low  as  possible 
so  the  value  of  AK  would  be  close  to  the  value  of  K.  A  minimum 
stress  of  .9  ksi  was  chosen  because  as  the  minimum  stress 
approached  zero,  the  load  wave  form  was  deformed  and  noise  or 
chatter  resulted  because  of  loose  fittings  and  slackness  in  the 
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machine.  To  prevent  slow  yielding  at  the  crack  tip  due  to  the  static 
load,  the  cyclic  load  was  not  stopped  during  the  crack  propagation 
stage  of  the  tests. 

3.4  Strain  Gauge  Instrumentation 

Strain  gauges  were  connected  in  such  a  fashion  that  only  the 
change  in  strain  due  to  the  presence  of  a  crack  was  recorded.  Five 
inch,  single  wire  strain  gauges  of  120  ohms  (compliance  gauges), 
were  placed  in  the  center  of  the  specimens  which  measured  the  average 
longitudinal  strain  in  the  cracked  section.  A  compensating  gauge 
of  a  quarter  inch  gauge  length  and  120  ohms  was  placed  near  the 
end  of  the  specimen  to  measure  the  average  longitudinal  strain  in 
the  uncracked  section  of  the  specimen. 

The  compensating  and  compliance  strain  gauges  were  placed 
on  each  side  of  the  sheet  specimen,  so  that  an  average  value  of  the 
strain  measurement  was  obtained.  The  compensating  gauges  were 
placed  in  opposite  arms  of  a  Wheatstone  bridge  circuit  while  the 
compliance  gauges  were  in  the  same  circuit  such  that  only  the  increase 
in  strain  due  to  the  presence  of  the  crack  was  measured.  Thus  the 
output  of  the  Wheatstone  bridge  circuit  was  twice  the  increase  of 
strain  due  to  the  crack.  A  variation  in  the  average  longitudinal 
stress  did  not  effect  the  compliance  measurement  because  the 
average  longitudinal  strain  was  compensated  for  in  the  Wheatstone 
bridge  circuit. 

Srawley,  Jones,  and  Gross  [15]  suggest  that  the  gauge  length 
should  be  four  times  the  length  of  the  longest  crack  in  the  SEN 
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specimen.  Since  the  gauge  length  of  five  inches  is  not  four  times 
the  longest  crack  considered  (1.55  inches),  a  slight  error  may  have 
resulted.  Errors  in  positioning  the  long  strain  gauges  in  the 
center  of  the  specimen  should  be  insignificant.  The  notch,  being 
only  on  one  side  of  the  SEN  specimen,  and  pin  loadings  caused  the 
SEN  specimen  to  bend  under  load.  Consequently  the  long  strain 
gauges  will  be  slightly  curved.  Srawley,  Jones,  and  Gross  [15] 
measured  the  bending  and  found  it  to  be  "small".  Thus  no  significant 
error  should  result  in  compliance  measurements  due  to  this  slight 
curvature  of  the  long  strain  gauges. 

An  SR-4  Strain  Indicator  and  an  oscilloscope  with  a  "Q"  unit 
was  used  to  analyze  the  output  of  the  strain  gauges.  The  SR-4 
Strain  Indicator  measured  the  mean  value  of  the  dynamic  strain  at 
the  center  of  the  specimen  due  to  the  crack.  The  accuracy  of  the 
SR-4  Strain  Indicator  was  adequate  for  compliance  calculations, 
whereas  the  accuracy  of  the  oscilloscope  was  not.  The  oscilloscope 
was  used  to  show  the  dynamic  strain  waveform  as  a  check  for  any 
distortion  or  resonance  of  the  strain  waveform.  This  was  necessary 
because  in  some  preliminary  tests,  the  dynamic  strain  waveform 
became  distorted  and  an  increase  in  magnitude  of  the  dynamic  strain 
of  up  to  25  percent  occurred  for  a  SEN  specimen  of  dimensions  0.05 
x  22  x  4  inches  over  a  relative  crack  length  of  0.35  to  0.415.  In 
the  tests  used  for  compliance  measurements,  only  the  0.125  inch  thick 
specimen  showed  a  similar  distortion.  The  distortion  was  noted  only 
for  a  short  interval  of  crack  growth,  thus  there  should  be  not  sig¬ 
nificant  effect  on  crack  propagation  results. 
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3.5  Crack  Length  ’leasurement 

A  measuring  microscope  was  used  to  both  observe  the  crack 
tip,  and  to  measure  the  crack  length  to  a  precision  of  a  thousandth 
of  an  inch.  A  30  power  microscope  was  used  in  initial  tests  and 
a  40  power  microscope  used  in  later  tests  made  the  crack  tip 
easier  to  define.  The  crack  tip  could  not  always  be  seen,  thus 
there  was  a  possible  human  error  present  in  defining  the  location 
of  the  crack  tip.  Erdogan  and  Roberts  [20]  found  the  sensitivity 
of  a  50  power  microscope  was  higher  than  the  accuracy  with  which 
the  crack  tip  could  be  defined,  as  some  contraction  in  thickness 
at  the  crack  tip  was  observed.  Thus  a  higher  power  microscope 
may  have  eliminated  errors  in  defining  the  crack  tip. 

A  stroboscope  allowed  the  crack  to  be  observed  in  the  fully 
opened  or  closed  states.  The  crack  tip  could  be  observed  best 
when  the  surface  of  the  specimen  was  left  in  the  oxidized  state 
with  a  minimum  of  scratches.  When  the  crack  was  in  the  fully 
opened  state,  the  flashing  light  made  the  shiny  surface  of  the 
crack  easy  to  see  against  the  darker  oxidized  surface  of  the  specimen. 

In  several  of  the  preliminary  tests,  the  alternating  load 
was  stopped  and  india  ink  injected  into  the  crack  near  the  crack 
tip.  The  ink  was  allowed  to  dry  in  the  crack  before  applying  the 
alternating  load.  Upon  complete  failure  of  the  specimens,  the 
shape  of  the  crack  front  could  be  seen  at  the  various  crack  lengths 
where  the  ink  had  been  injected  and  dried. 

Observations  from  these  preliminary  tests  were  used  to  calculate 
the  crack  lengths  at  the  center  portion  of  the  specimen. 
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CHAPTER  IV 


TEST  RESULTS  AND  DISCUSSION 


4.1  Compliance  and  Stress  Intensity  Factor  Curves 

The  original  test  data,  that  is  the  compliance  measurements 
and  surface  crack  lengths,  is  presented  in  tabular  and  graphical 
form.  The  stress  intensity  factors  are  presented  in  graphical 
form  for  comparison.  The  differences  of  the  results  are  discussed 
with  respect  to  different  characteristics  of  the  tests  and  possible 
sources  of  error. 

The  set  of  crack  lengths,  compliance  readings,  and  number  of 
loading  cycles  are  recorded  for  six  tests  in  Tables  II  to  VII. 

The  first  five  crack  lengths  are  sawcuts  of  approximately  0.006  inches 
in  width  and  the  remaining  are  fatigue  crack  lengths.  The  compliance 
readings  for  the  sawcuts  provide  points  on  the  curves  below  a 
relative  crack  length  of  0.2.  The  useful  region  for  fatigue  testing 
on  the  SEN  specimen  is  between  0.2  and  0.5  relative  crack  lengths. 

The  compliance  readings  and  relative  crack  lengths  are  plotted  for 
the  tests  in  Figures  4.1,  4.2,  and  4.3.  To  avoid  confusion,  only 
about  one  half  of  the  readings  given  in  the  tables  are  shown  in  the 
figures . 

The  compliance  readings  for  several  tests  of  Srawley,  Jones, 
and  Gross  [15]  are  shown  in  Figure  4.5.  They  used  the  average 
compliance  readings  for  two,  twelve  inch  long  specimens  to  calculate 
G  and  to  be  used  for  publication  in  later  papers  [24].  The  average 
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compliance  readings  of  these  two  specimens  are  used  for  comoarison 
in  Figures  4.1,  4.2,  and  4.3.  The  difference  between  tests  and  the 
scatter  in  each  test  are  probably  the  result  of  errors  in  crack 
length  measurements,  errors  in  compliance  readings,  and  the 
fatigue  crack  changing  from  a  flat  to  a  slanted  crack. 

Compliance  readings  for  several  tests  had  to  be  disregarded 
because  of  errors  due  to  the  binding  of  specimens  in  the  holding 
blocks  and  large  temperature  difference  during  tests.  The  accuracy 
of  the  travelling  microscope  was  adequate  for  accurate  crack 
length  measurements.  A  small  error  may  be  present  in  crack  length 
measurements  due  to  the  difficulty  of  seeing  the  exact  location 
of  the  crack  tip,  but  would  not  exceed  one  hundredth  of  an  inch. 

From  the  tabulated  data,  it  can  be  seen  that  the  load  was  applied 
a  few  thousand  times  before  a  notch  was  cut.  This  initial  loading 
served  as  a  shakedown  period,  required  to  get  stable  compliance 
readings  for  a  zero  crack  length.  All  the  compliance  readings 
were  taken  when  the  specimen  was  under  the  dynamic  load  as  there 
was  an  unexplainable  difference  between  the  compliance  reading  of 
the  mean  dynamic  load  and  the  compliance  reading  of  the  static 
load.  The  difference  was  fairly  constant  for  the  entire  test. 

A  total  of  twenty-one  SEN  specimens  were  tested.  Most  of  the 
preliminary  tests  had  large  errors  in  the  compliance  readings.  Changes 
in  the  experimental  methods  such  as;  use  of  compensating  gauges, 
use  of  a  stroboscope,  control  of  temperature,  and  accurate  setting 
of  end  fixtures  eliminated  most  of  the  errors.  As  seen  from  the 
scatter  of  points  and  differences  for  tests  on  identical  specimens, 
there  must  still  be  small  errors. 
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The  least  squared  error  method  was  used  to  fit  the  compliance 
versus  relative  crack  length  curves  to  polynomial  expressions  as 
discussed  in  Section  2.7.  The  stress  intensity  factor  curves 
obtained  from  the  compliance  curves  are  shown  in  Figures  4.5  to  4.10 
for  specimens  of  thickness  0.050,  0.080,  and  0.125  inches.  The 
stress  intensity  factor  curves  for  the  plane  stress  problem  as 
calculated  from  Srawley,  Jones,  and  Gross'  [15]  work  is  shown  on 
each  figure  for  comparison.  Their  K  curve  is  also  shown  for  the 
plane  strain  condition,  as  well  as  the  plane  stress  condition  in 
Figure  4.11.  Comparing  the  curves  over  a  relative  crack  length 
0.2  to  0.45,  tests  numbered  1,  3,  4,  and  6  are  slightly  above 
Srawley,  Jones,  and  Gross'  curves;  whereas  tests  numbered  2  and  5 
are  slightly  below  those  of  Srawley,  Jones,  and  Gross.  From  these 
observations,  it  appears  that  the  K  curve  of  Srawley,  Jones,  and 
Gross  is  applicable  to  fatigue  cracks  in  SEN  specimens  if  the 
plastic  zone  correction  factor  is  applied  to  the  fatigue  crack 
length. 

An  example  of  the  effect  of  the  plastic  zone  correction  on 
K  can  be  obtained  from  Figure  4.12  which  shows  K  calculated  with 
and  without  the  plastic  zone  and  crack  curvature  corrections.  There 
is  a  significant  shift  of  the  K  curve  without  the  use  of  this 
correction.  There  are  no  notable  differences  in  the  curves  for 
specimens  of  different  thicknesses.  In  all  test  curves,  except 
test  number  4,  the  value  of  the  stress  intensity  factor  increased 
slightly  in  relation  to  that  of  Srawley,  Jones,  and  Gross,  as  a 
relative  crack  length  of  0.45  was  approached.  The  slant  of  the 
crack  surface  increased  as  the  relative  crack  length  increased. 


36 


This  is  ullustrated  in  Figure  4.13.  The  slanting  crack  surface  is 
the  only  observed  factor  that  could  possibly  cause  the  increased 
value  of  K  relative  to  the  value  of  K  as  obtained  by  Srawley,  Jones 
and  Gross.  This  increase  would  not  amount  to  five  percent,  thus 
it  is  not  too  significant. 

4.2  Transition  of  Crack  Surfaces 

The  cross-sections  of  test  specimens  are  shown  in  Figure  4.13 
for  various  relative  crack  lengths.  The  fatigue  crack  surface  is 
initially  flat  at  the  start  of  the  fatigue  crack  and  slanted  near 
the  center  of  the  specimen.  The  surface  changes  by  the  formation 
of  a  "shear  lip"  growing  from  the  surface  of  the  specimen  across 
the  crack  to  the  center,  at  which  time  the  crack  surface  would  be 
at  a  fairly  constant  slope  of  about  20  degrees.  In  earlier  tests, 
the  stress  ratio  (R)  was  larger  and  a  slope  of  about  38  degrees 
was  obtained.  As  seen  from  Figure  4.13  the  slope  does  not  change  uni¬ 
formly.  In  the  thinner  specimens  (0.050  inches  thick),  the  slanted 
surface  developed  rather  "fast"  and  shifted  back  and  forth.  In  the 
thicker  specimens  (0.125  inches  thick),  the  slanted  surface  developed 
slowly  and  more  uniformly. 

At  times  a  "V"  crack  surface  developed  rather  than  a  slanted 
surface.  The  "V"  crack  caused  the  crack  to  deviate  from  the  straight 
crack  path  in  a  direction  opposite  to  the  point  of  the  "V". 

Irwin  [21]  has  studied  the  transition  of  the  crack  surface  for 
fracture  toughness  tests.  He  found  from  equations  3a,  3b,  and  3c 
that  the  distance  r  from  the  crack  tip  at  which  the  predicted  tensile 
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stress  equals  the  yield  stress,  is  greatest  at  0  =  60  degrees  where 

0 . 84  K2 


r  = 


TT 


(23) 


ys 


When  the  plate  thickness  B,  was  equal  to  twice  r  of  equation  23, 
Irwin  concluded  that  the  transition  of  the  crack  from  flat  to  slant 
was  beginning,,  That  is 


B  =  2r  and  (24) 

/Btt  a 

K  =  - (25) 

/1 .68 


He  showed  that  there  was  good  agreement  for  the  results  of  fracture 
toughness  tests.  Applying  equations  24,  25,  and  K  curves  to  test 
specimens,  the  relative  crack  lengths  at  which  the  surface 
transition  would  occur  should  be  0.35,  0.38,  and  0.45  for  specimens 
of  thickness  0.050,  0„080,  and  0.125  inches  respectively.  As  seen 
in  Figure  4.13  the  transition  occurs  at  smaller  relative  crack 
lengths . 

Hertzberg  and  Paris  [19]  use  the  same  reasoning  as  Irwin  [21] 
for  transition  of  fatigue  cracks  on  thinner  specimen  except  that 


K2 
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and 

(26) 

/ttF  0 

ys 

(27) 

They  found  the  K  of  equation  27  occurred  approximately  at  the  point 
when  the  flat  fracture  was  completely  eliminated.  Using  this 
reasoning,  the  transition  should  occur  at  relative  crack  lengths  of 
0.28,  0.32,  and  0.41  for  0.050,  0.080,  and  0.125  inch  thick  specimens, 
respectively.  These  crack  lengths  agree  better  than  those  obtained 
from  equation  25  for  the  cross  section  on  Figure  4.13, 
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4.3  Validity  of  the  Stress  Intensity  Factor 

The  stress  intensity  factor  was  obtained  from  the  strain 
energy  release  rate,  This  rate  is  developed  from  energy  concepts 
for  cracks  with  no  restrictions  as  to  a  perfectly  elastic  or 
plastic  solid,  or  a  slanted  or  flat  crack.  The  stress  intensity 
factor  has  been  developed  for  a  flat  crack  in  a  perfectly  elastic 
solid  under  conditions  of  plane  stress  or  plane  strain,,  The  G 
expressions  obtained  for  the  test  specimens  are  completely  correct 
according  to  the  definition  of  G„  The  K  expressions  obtained  for 
the  fatigue  cracks  are  not  correct  according  to  the  assumptions 
for  which  K  was  defined.  The  effects  of  the  plastic  zone  has  been 
approximately  corrected  for  when  calculating  K  but  no  correction 
could  be  applied  for  the  slanting  of  the  crack  tip. 

A  slant  crack  tip  is  under  a  shear  and  tension  loading.  Stress 
intensity  factors  have  been  developed  for  shear  and  tension  loads 
on  flat  crack  tips  and  could  possibly  be  related  to  slanted  cracks. 
Paris  and  Si h  [14]  have  related  G  to  K  for  a  shear  load  in  a  manner 
similar  to  that  relating  G  to  K  for  a  tension  load  as  in  Section  2.4. 
Since  only  one  value  of  G  was  found  experimentally,  there  is  no 
way  of  knowing  what  portion  of  G  would  be  due  to  the  shear  loading 
or  to  the  tension  loading  for  a  slanted  crack0  The  K  defined  for 
tension  loadings  has  been  accepted  as  a  parameter  that  is  correlated 
with  the  crack  propagation  rate  [19,  20],  but  it  should  be  kept  in 
mind  that  this  K  is  not  an  accurate  description  of  the  stress  fields 
for  a  slanted  crack  but  may  occasionally  be  encountered. 
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TABLE  II 

Experimental  Data  for  a  0.050  Inch  SEN  Specimen 

Test  No.  1 


Number  of 
load  cycles 

1 000 1 s 

Crack 
length 
i  n , 

EAC 

20 

0 

0 

21 

.12 

.016 

22 

.23 

.030 

22 

.37 

.056 

23 

.50 

.096 

24 

,635 

.142 

72 

.648 

.153 

80 

.656 

.158 

90 

.665 

.156 

100 

.679 

.172 

110 

.694 

.178 

120 

.704 

.184 

130 

.715 

.195 

140 

,734 

.200 

150 

,752 

.219 

160 

765 

.233 

170 

,788 

.249 

180 

.800 

.263 

190 

.826 

.285 

200 

.848 

.298 

210 

,865 

.318 

220 

,882 

.353 

230 

.909 

.372 

240 

,940 

.400 

245 

.956 

.425 

Number  of  Crack 

loading  cycles  length 
1000's  in. 


250 

.969 

255 

.987 

260 

1.003 

265 

1,032 

270 

1.048 

275 

1.075 

280 

1.095 

285 

1.118 

290 

1  .143 

295 

1  .179 

300 

1.204 

305 

1  ,238 

309 

1  ,259 

313 

1,301 

316 

1  328 

319 

1,355 

322 

1,385 

325 

1  ,424 

328 

1  .465 

331 

1.505 

334 

1,560 

336 

1  .610 

338 

1.64 

339 

1  ,67 

341 

1,71 

342 

1.75 

343 

1  .85 

344 

Fai led 

EAC 


.446 

.462 

.487 

.515 

.554 

.583 

.616 

.656 

.715 

.783 

.819 

.916 

.990 

1.066 

1.150 

1  .236 
1  .335 
1.460 
1.610 
1.811 

2.101 

2.319 

2.801 
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TABLE  III 

Experimental  Data  for  a  0.050  Inch  SEN  Specimen 

Test  No.  2 


Number  of 
load  cycles 
1000's 

Crack 
length 
i  n  o 

Number  of 

EAC  loading  cycles 

1 000 1 s 

Crack 

length 

in. 

EAC 

50 

0 

0 

235 

.946 

.394 

50 

.12 

.016 

245 

.988 

.443 

50 

.23 

.019 

250 

1 .003 

.475 

50 

.37 

,042 

255 

1.032 

.496 

50 

.49 

.081 

260 

1 .067 

.550 

56 

.617 

.136 

265 

1 .082 

.575 

95 

.627 

.147 

270 

1.112 

.609 

105 

.638 

.147 

275 

1.140 

.660 

115 

.654 

.151 

280 

1 .182 

.745 

125 

.669 

.159 

284 

1.213 

.785 

135 

.686 

.173 

288 

1  .250 

.853 

145 

.702 

.189 

292 

1  .283 

.956 

155 

.722 

.202 

296 

1.327 

1 .160 

165 

.737 

.205 

300 

1.373 

1  .222 

175 

.765 

.218 

303 

1.411 

1.346 

185 

.794 

.250 

306 

1.462 

1  .50 

195 

.815 

.,269 

308 

1.494 

1.63 

205 

.844 

.296 

310 

1.542 

1.80 

215 

.880 

.  322 

312 

1  .587 

2.05 

225 

.908 

.338 

313 

1  .612 

2.20 

314 

1.629 

2.30 
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TABLE  IV 

Experimental  Data  for  a  0.080  Inch  SEN  Specimen 

Test  No.  3 


Number  of 

Crack 

Number  of 

Crack 

load  cycles 

length 

EAC 

loading  cycles 

length 

EAC 

1000's 

in. 

1000*  s 

in. 

33 

.00 

0 

240 

,999 

.512 

33 

o  09 

.01 3 

245 

1.020 

.537 

37 

o  19 

.024 

250 

1.030 

.562 

41 

.31 

.046 

255 

1.048 

.594 

46 

.40 

.075 

260 

1.068 

.620 

50 

51 

.101 

265 

1.097 

.675 

54 

.59 

.148 

270 

1.138 

.742 

90 

.623 

o  1 56 

2.75 

1.165 

.798 

110 

,653 

.181 

280 

1.196 

.865 

124 

.675 

CO 

o 

CVI 

o 

285 

1.241 

.976 

135 

o  701 

.214 

287 

1.278 

1 .100 

140 

o  71 3 

.220 

291 

1 .316 

1.182 

150 

o  725 

.238 

294 

1  .351 

1.278 

160 

.745 

.249 

297 

1  .379 

1.390 

170 

.761 

.268 

299 

1.409 

1.521 

180 

.790 

.290 

300 

1 .429 

1.588 

190 

.820 

.312 

302 

1.454 

1.675 

200 

.847 

.353 

303 

1  471 

1.791 

205 

.872 

.364 

306 

1 .55 

210 

.888 

.386 

314 

1.9 

215 

.903 

,404 

220 

.925 

,428 

225 

,937 

441 

230 

.959 

.462 

235 

.974 

.485 

I  ' 


■ 

. 
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TABLE  V 

Experimental  Data  for  a  0.080  Inch  SEN  Specimen 

Test  No.  4 


Number  of 
load  cycles 
1000's 

Crack 

length 

inn 

EAC 

Number  of 
loading  cycles 
1000's 

Crack 

length 

in. 

EAC 

19 

0 

0 

245 

1,013 

.530 

20 

.1 

.013 

250 

1.039 

.567 

21 

.25 

.033 

255 

1.067 

.655 

22 

o375 

.063 

260 

1.105 

.675 

23 

.49 

.098 

265 

1.138 

.739 

26 

,605 

.140 

270 

1.170 

.835 

70 

.610 

.142 

274 

1.205 

.904 

80 

o616 

.151 

278 

1.228 

.981 

90 

,625 

.156 

282 

1.290 

1.110 

100 

.631 

.162 

285 

1.329 

1.211 

111 

„  650 

.169 

288 

1.375 

1.328 

120 

.657 

.180 

291 

1,415 

1.491 

130 

.678 

.202 

294 

1  472 

1.710 

140 

.695 

.208 

297 

1 . 526 

1.930 

150 

.716 

.219 

300 

1.620 

2.161 

160 

,734 

.230 

304 

175 

170 

758 

.249 

305 

1 . 85 

180 

780 

.271 

190 

.805 

.293 

200 

.828 

.322 

210 

.863 

.348 

220 

.889 

,383 

230 

,934 

.427 

235 

.959 

,461 

240 

993 

,497 
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TABLE  VI 

Experimental  Data  for  a  0.125  Inch  SEN  Specimen 

Test  No.  5 


Number  of 
load  cycles 
1000's 

Crack 

length 

in. 

Number  of 

EAC  loading  cycles 
1000 ' s 

Crack 

length 

in. 

EAC 

35 

0 

0 

206 

1.141 

.659 

40 

.012 

.016 

209 

1 .185 

.714 

45 

.022 

.024 

212 

1.220 

.765 

50 

.036 

.057 

214 

1.249 

.816 

55 

.050 

.093 

216 

1.290 

.912 

56 

.062 

.128 

219 

1.339 

1.072 

99 

.65 

.131 

222 

1.384 

1.190 

109 

.676 

.145 

223 

1.411 

1.285 

119 

.693 

.166 

225 

1.464 

1.462 

129 

.721 

.184 

227 

1.511 

1.685 

139 

.753 

.205 

228 

1.554 

149 

.784 

.230 

229 

1.602 

159 

.813 

.255 

230 

1.65 

169 

.874 

.301 

230.5 

1.75 

179 

.927 

.348 

231 

1.90 

184 

.945 

.386 

189 

.984 

.427 

194 

1.024 

.473 

200 

1.088 

.576 

203 

1.106 

.615 
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TABLE  VII 

Experimental  Data  for  a  0.125  Inch  SEN  Specimen 

Test  No.  6 


Number  of 
load  cycles 
1000 1 s 

Crack 
length 
in . 

Number  of 

EAC  loading  cycles 
1000 1 s 

Crack 

length 

in. 

EAC 

49 

0 

0 

185 

.962 

.428 

50 

.12 

.011 

190 

.994 

.468 

51 

.24 

.030 

195 

1.027 

.523 

52 

.38 

.065 

200 

1.067 

.595 

56 

.612 

.156 

203 

1.102 

.635 

62 

.612 

.136 

206 

1.129 

.687 

98 

.623 

.145 

209 

1.148 

.741 

108 

.669 

.156 

212 

1.190 

.823 

119 

.675 

.181 

215 

1.225 

.872 

129 

.706 

.194 

218 

1.272 

1.015 

139 

.740 

.206 

220 

1 .300 

1.090 

149 

.777 

.240 

222 

1.335 

1.190 

159 

.815 

.277 

224 

1  .384 

1.380 

170 

.869 

.324 

226 

1.435 

1.582 

180 

.926 

.383 
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Figure  4.1  Compliance  Curve  for  0.050  inch  Thick  SEN  Specimen 
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Figure  4.2  Compliance  Curves  for  0.080  Inch  SEN  Specimen 
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Figure  4.3  Compliance  Curves  for  0.125  Inch  SEN  Specimen 
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Figure  4.4  Compliance  Curves  of  Srawley,  Jones,  and  Gross  [15] 
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Figure  4.5  Dimensionless  Stress  Intensity  Factor  Plot  for  Test  #1 


DIMENSIONLESS  STRESS  INTENSITY  FACTOR  -  rC/a/T 
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Figure  4.6  Dimensionless  Stress  Intensity  Factor  Plot  for  Test  #2 


DIMENSIONLESS  STRESS  INTENSITY  FACTOR  -  K/a/w 
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Figure  4.7  Dimensionless  Stress  Intensity  Factor  Plot  for  Test  #3 
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Figure  4.8  Dimensionless  Stress  Intensity  Factor  Plot  for  Test  #4 
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Figure  4.9  Dimensionless  Stress  Intensity  Factor  Plot  for  Test  #5 


DIMENSIONLESS  STRESS  INTENSITY  FACTOR  -  K/o/B" 
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Figure  4.10  Dimensionless  Stress  Intensity  Factor  Plot  for  Test  #6 


DIMENSIONLESS  STRESS  INTENSITY  FACTOR  -  K/o/vT 
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Figure  4.11  Dimensionless  Stress  Intensity  Factor  Plot  for 

Srawley,  Jones,  and  Gross  [15]  12  Inch  Long  Specimen 


DIMENSIONLESS  STRESS  INTENSITY  FACTOR  -  K/oJTT 
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Figure  4.12  Dimensionless  Stress  Intensity  Factor  Plots  With 

and  Without  the  Plastic  Zone  Corrections  for  Test  #6 


0.050  IN.  THICK  SPECIMENS 


57 


Test  #1 

0 


Test  #2 


0.080  IN.  THICK  SPECIMENS 


Test  #3 


Test  #4 


0.126  IN.  THICK  SPECIMENS 


Test  #5 


RELATIVE  CRACK  LENGTH  a/w 
POSITION  OF  CROSS-SECTIONS 


Figure  4.13  Cross-Sections  of  Crack  Surfaces 
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4.4  Stress  Intensity  Factor  Correlation  With  the  Logarithm  of 
the  Crack  Propagation  Rate. 

The  fatigue  crack  propagation  rate  has  been  considered  to  be 
a  function  of  AK  and  several  types  of  plots  have  been  used  for  the 
correlation.  Since  there  is  a  great  deal  of  scatter  using  linear 
coordinate  plots  of  AK  versus  da/dN,  logarithm  coordinates  have  been 
used.  The  AK  versus  log  (da/dN)  plot  has  been  used  for  close 
observation  of  test  data  while  the  log  (AK)  versus  log  (da/dN)  plot 
has  been  used  to  show  general  agreement  of  test  data  from  various 
sources  for  a  specific  material. 

The  AK  versus  log  (da/dN)  plot  for  aluminum  alloys  can  be 
described  as  having  three  branches:  a  lower,  central,  and  upper 
branch  as  illustrated  in  Figures  4.14  to  4.18.  The  lower  branch 
of  the  curve  occurs  when  AK  is  below  a  value  of  6000  Ib(in)1^  as 
illustrated  in  Figures  4.15  and  4.17.  The  data  for  the  lower  part 
of  the  curve  indicates  a  fluctuation  in  the  crack  propagation  rate 
as  the  crack  grows.  This  is  possibly  due  to  the  crack  growth  rate 
being  discontinuous  and/or  the  accuracy  of  crack  length  measurements 
not  being  sufficient  for  an  accurate  determination  of  da/dN.  A 
line  through  the  data  points  of  the  lower  branch  has  a  very  gentle 
slope.  The  crack  surfaces  are  usually  completely  flat  for  this 
stage  of  the  crack  growth. 

The  central  branch  of  the  curve  is  indicated  in  Figures  4.14 
to  4.18.  A  straight  line  can  be  drawn  through  the  data  points  for 
this  branch  of  the  curve.  A  great  deal  of  scattering  can  be  seen 
in  the  curves,  and  is  even  more  pronounced  when  the  stress  ratio  R 
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is  of  a  lower  value,  as  for  example  in  Figure  4.14  and  4.15.  The 
slope  and  position  of  the  central  branch  of  the  curve  depends  on 
the  specimen  thickness.  The  thinner  specimens  have  slower  propagation 
rate  than  the  thicker  specimens  for  the  same  value  of  AK  as  seen 
from  the  above  figures.  Increasing  the  stress  ratio  shifts  the 
curve  to  the  right  as  seen  in  Figures  4.14  and  4.15  and  thus  increases 
the  crack  propagation  rate.  The  crack  surface  is  in  a  transition  phase 
changing  from  flat  to  slanted  over  the  central  branch  of  the  curve. 

The  upper  branch  of  the  correlation  is  shown  in  Figure  4.16 
and  begins  at  a  point  where  the  central  portion  deviates  to  the 
right.  The  location  of  this  upper  portion  depends  upon  the  stress 
ratio.  Increasing  the  stress  ratio  will  cause  this  portion  of  the 
curve  to  begin  at  a  lower  value  of  AK„  Thus  this  portion  probably 
depends  more  on  K  than  AK. 

The  AK  versus  log  (da/dN)  plot  can  be  used  to  illustrate  the 
effects  of  test  variables  quite  well.  In  the  figures,  the  effects 
of  specimen  thicknesses,  stress  ratios,  and  annealing  on  the  crack 
propagation  ratio  can  be  seen. 

4.5  Relating  the  Correlation  of  the  Stress  Intensity  Factor  and 

Crack  Propagation  Rate  to  Electron  Fractography 

Piper,  Quist,  and  Anderson  [25]  have  plotted  data  for  the 
aluminum  alloy  7178-T6,  using  the  AK  versus  log  (da/dN)  plot  and 
they  discuss  the  correlation  in  three  branches  as  in  Section  4.4. 
Electron  mi crofractography  studies  of  the  crack  surfaces  was  used 
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Figure  4.14  Crack  Propagation  Curve  for  Test  #1  to  #6  on  SEN  Specimen 
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Figure  4.15  Crack  Propagation  Curve  for  Preliminary  Test  on  SEN  Specimen 
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Figure  4.16  Crack  Propagation  Curve  Showing  "Upper 


crniax  =  3.5  k.s.i.  R=.5  Freq.  =  1800  cpm  Width 
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Figure  4.18  Crack  Propagation  Curve  for  Preliminary  Tests  on  Center-Notch  Specimen 
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to  discuss  the  "mechanism"  of  crack  growth  for  each  portion  of  the 
curve.  For  the  lower  portion,  fatigue  striations  were  evident  over 
all  the  replicated  area  of  the  crack  surface.  They  were  also  evident 
over  most  of  the  replicated  area  of  the  crack  surface  of  the 
central  portion  and  only  a  small  percentage  of  the  replicated  area 
of  the  crack  surface  of  the  upper  portions.  They  describe  the 
remaining  replicated  area  as  represented  by  internal  voids  which 
coalesce  by  a  process  of  internal  and  localized  necking. 

Christensen  and  Harmon  [12]  compare  the  striations  to  cyclic 
loadings  and  conclude  that  at  higher  values  of  K  the  crack 
propagates  by  means  other  than  producing  striations. 

Hertzberg  and  Paris  [19]  studied  the  striations  on  fatigue 
surfaces  of  2024-T3  specimens  of  thicknesses  0.063,  0.094,  and  0.126 
inches.  They  concluded  that  fatigue  striations  were  present  in 
plane  strain  regions  and  absent  from  plane  stress  regions.  Striations 
completely  covered  the  flat  crack  area  and  only  part  of  the  area  of 
the  slanted  crack.  They  described  the  remaining  area  of  slanted 
cracks  as  having  examples  of  ductile  ruptures  called  elongated 
dimples,  which  have  a  long  axis  perpendicular  to  the  direction  of 
crack  propagation.  The  striation  agreed  with  the  loading  cycles 
for  a  crack  propagation  rate  of  1  to  100  microinches  per  cycle 
which  corresponds  to  the  central  portion  of  the  AK  versus  log  (da/dN). 

From  these  observations,  it  appears  that  the  crack  propagates 
in  a  manner  producing  striations  for  the  lower  and  central  portions 
of  the  AK,  log  (da/dN)  curve.  For  the  lower  portion  the  cyclic  load 
probably  does  not  produce  striations  continuously  while  in  the  upper 
portion  of  the  curve  the  crack  probably  propagates  by  a  mechanism 
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other  than  that  which  produces  striations.  Tetelman  and  McEvily  [10] 
state  that  any  model  used  for  fatigue  crack  propagation  must  account 
for  the  striations,  and  other  surface  marks,  and  the  observation 
of  plastic  deformation  taking  place  at  the  crack  tip. 


4.6  Logarithm  of  Stress  Intensity  Factor  Correlated  with  Logarithm 
of  Crack  Propagation  Rate 


Paris  and  Erdogan  [18]  have  used  the  log  (AK)  versus  log  (da/dN) 
plot  to  show  agreement  for  crack  propagation  data  of  2024-T3  material 
from  various  sources.  Data  for  crack  propagation  rates  of  less  than 
one  micro  inch  per  cycle  to  10,000  micro  inches  per  cycle  was 
plotted  on  the  one  graph.  A  straight  line  of  slope  1/4  appeared 
to  fit  data  for  the  wide  range  of  crack  propagation  rates.  The 
equation  of  the  line  being 

m  -  (28) 


where  M  is  a  material  constant.  Tentatively,  equation  28  could 
be  used  to  predict  crack  propagation  but  it  basically  shows  that 
larger  changes  in  the  crack  propagation  rate  may  occur  for  relatively 
small  changes  in  AK. 

More  recently  Roberts  and  Erdogan  [20]  have  presented  a 
crack  propagation  equation  which  accounts  for  the  effect  of  the  stress 
ratio  R  as  well  as  AK  on  the  crack  propagation  rates  for  fluctuating 
tensile  loading  and  bending  loadings.  They  initially  assumed  that 
da/dN  should  depend  on  the  plastic  deformation.  Relating  the  plastic 
deformation  to  the  stress  intensity  factors,  they  found  that  the 
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crack  propagation  rate  must  be  dependent  on  the  stress  intensity 
factors.  They  assumed  the  crack  propagation  equation  to  be: 
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(29) 


where  B  is  a  material  constant.  From  da/dN  tests  on  2024-T3  and 
7075-T6  materials,  they  concluded  the  general  da/dN  law  to  be: 


$]•  “  B(K)2  (AK)2 

w  -  B<^>2  (AK)4 


(30) 


(31) 


or 


The  specific  equation  that  represents  the  test  data  of  2024-T3 
specimen  under  fluctuating  tension  is: 


^  =  2.68  x  10"19  (y^O 1  * 54  ( AK) 3  * 62 


(32) 


Data  of  tests  1,  4,  and  6  are  plotted  using  log  AK  versus  log  da/dN 
coordinate  in  Figure  4.19  with  the  equation  of  each  curve  shown  on 
the  figure.  The  lines  of  equations  31  and  32  are  shown  on  the  figure 
for  comparison.  The  test  data  curves  are  above  the  curves  of 
equation  32,  but  of  approximately  the  same  slope.  The  difference 
may  be  due  to  the  higher  frequency  of  the  tests  or  other  test 
parameters  that  affect  da/dN.  The  data  of  the  2024-T3  material  of 
Roberts  and  Erdogan  [20]  agrees  fairly  well  with  that  of  Hertzberg 
and  Paris  [19]  but  not  with  the  test  results  of  tests  1  to  6. 
Hertzberg  and  Paris  [19]  indicated  that  higher  loading  frequencies 
slowed  down  the  crack  propagation  rate.  Tests  numbered  1  to  6 
had  a  loading  frequency  of  1800  cpm  while  those  of  references  19 
and  20  were  at  140  cpm.  The  yield  stress  of  the  material  tested 
varied  between  52  to  57  k.s.i.  which  was  a  little  higher  than  the 
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yield  stress  of  the  material  used  by  Hertzberg  and  Paris  [19] 

(45  k.s.i.).  The  yield  stress  was  raised  to  the  fourth  power  and 
included  as  part  of  the  material  constant  in  Roberts  and  Erdogan's 
work.  The  difference  in  the  yield  stress  may  be  due  to  differences 
in  properties  of  the  material  in  manufacturing  and/or  to  age 
hardening.  The  age  of  the  2024-T3  material  is  unknown  and  its 
yield  stress  might  have  increased  due  to  the  age  hardening  process. 

According  to  Piper,  Quist,  and  Anderson  [25],  slight 
variations  of  chemical  composition  may  change  the  mechanical 
properties  more  than  variations  in  mill  practices.  They  showed 
variations  of  the  chemical  composition  of  the  7000  series  of 
aluminum  alloys,  greatly  affected  the  crack  propagation  rates 
and  fracture  toughness  values. 
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CHAPTER  V 


CONCLUSIONS 


5.1  Stress  Intensity  Factor 

The  stress  intensity  factor  which  was  found  using  a  flat 
fatigue  crack  for  a  notch  and  a  plastic  zone  correction,  agreed 
with  the  stress  intensity  factors  obtained  by  Srawley,  Jones  and 
Gross  [15],  and  Gross,  Srawley,  and  Brown  [14]  using  an  artificial 
ideal  elastic  crack.  The  K  found  for  a  slanted  crack  was  slightly 
greater  than  the  K  for  a  flat  crack.  The  application  of  K  to  a 
slanted  crack  is  not  completely  correct  since  K  was  defined  for 
a  flat  crack  under  tension  loading  at  the  crack  tip. 

The  thickness  of  the  specimens  did  not  affect  the  stress 
intensity  factor  for  specimens  of  thickness  0.05,  0.080,  and  0.125. 

The  use  of  long  strain  gauges  was  found  to  be  an  accurate 
method  of  measuring  the  SEN  specimen  compliance. 

5.2  Crack  Propagation 

The  crack  propagation  rate,  for  the  bare  aluminum  alloy 
2024-T3,  was  found  to  be  a  function  of  the  stress  intensity  factor, 
AK,  with  many  factors  affecting  this  function.  The  AK  versus 
log  (da/dN)  plot  was  used  to  show  the  effect  of  single  factors  on 
the  da/dN,  AK  function.  The  lower  part  of  the  plot  was  unaffected 
by  specimen  thickness  whereas  for  the  other  parts  of  the  plot, 
the  thinner  specimens  (0.050  inches)  had  a  slower  crack  propagation 
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rate  than  the  thicker  specimens  (0.125  inches)  for  the  same  value 
of  AK.  The  log  AK  versus  log  da/dN  plot  was  found  useful  in 
showing  agreement  of  a  wide  range  of  data  for  one  material  and  for 
comparison  of  crack  propagation  properties  of  various  materials. 

Increasing  the  stress  range  caused  increases  in  the  crack 
propagation  rate  for  a  given  value  of  AK.  Roberts  and  Erdogan  [20] 
incorporated  the  effects  of  R  in  the  crack  propagation  equation. 

The  crack  propagation  rate  was  written  as  a  function  of  R  and  AK. 
From  tests,  the  exponent  of  the  AK  term  in  the  da/dN  equation  was 
found  to  be  4  and  agreed  with  that  of  Roberts  and  Erdogan  [20]  and 
Hertzberg  and  Paris  [18]. 

The  2024-T3  aluminum  alloy  used  in  tests  had  a  slower  crack 
propagation  rate  than  that  of  references  20  and  18,  which  caused 
a  significant  difference  of  the  material  constants  in  crack 
propagation  equation.  Reasons  for  the  slower  crack  propagation  rate 
might  be  attributed  to  the  higher  frequency  of  the  cyclic  loading 
(1800  cpm  as  compared  to  140  cpm)  and  to  the  higher  yield  strength 
(56,000  psi  as  compared  to  45,000  psi). 

The  transition  of  a  fatigue  crack  surface  from  being  flat 
to  slanted  appears  to  depend  on  the  plastic  zone  size  and  is  an 
indication  of  a  change  from  plane  strain  condition  to  a  plane  stress 
condition.  The  transition  for  thinner  specimens  (0.050  inches) 
was  more  irregular  than  that  of  thicker  specimens  (0.125  inches) 
and  occurred  over  a  shorter  period  of  crack  growth. 

The  stroboscope  and  high  power  microscope  were  successfully 
used  to  enable  continuous  measurement  of  the  crack  length  during 
testing. 
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